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Abstract 

We consider equations of the form 

(D(p)n) (t) = -Xu{t), t > 0, 

where A > 0, D(p) is a distributed order derivative, that is 

1 



D^") is the Caputo-Dzhrbashyan fractional derivative of order a, p is a positive measure. 

The above equation is used for modeling anomalous, non-exponential relaxation pro- 
cesses. In this work we study asymptotic behavior of solutions of the above equation, 
depending on properties of the measure p. 
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1 INTRODUCTION 



Anomalous, non-exponential, relaxation processes occur in various branches of physics; see e.g. 
[SlIIHlEI] and references therein. Just as the exponential function e"'^* (A > 0) appearing in the 
description of classical relaxation is a solution of the simplest differential equation u' = —Xu 
(with the initial condition u{0) = 1), new kinds of derivatives are used to obtain models of slow 
relaxation. 

It is generally accepted now that the power law of relaxation corresponds to the Cauchy 
problem 

(t) = -Xu{t), t > 0; u{0) = 1, (1) 

where 

d 



Til -a) 



dt 



{t-Tyu{T) dr-r^uiO) 



(2) 



is the Caputo-Dzhrbashyan fractional derivative of order a G (0, 1); we refer to [7l[l3] for various 
notions and results regarding fractional differential equations. The solution of the problem (1) 
has the form u{t) = £'„(— At") where is the Mittag-LefHer function, and u{t) ~ Ct~" (here 
and below we denote various positive constants by the same letter C), as t — > oo. This kind 
of evolution describes the temporal behavior related to a-fractional diffusion typical for fractal 
media. 

A still slower, logarithmic relaxation [TTl [TB] is described by the Cauchy problem 

(D('^)m) (t) = -Xu{t), t > 0; m(0) = 1, (3) 

where 

1 

{B^^^u) (t) = [ (t)/i(a) da, (4) 



/i is a non-negative continuous function on [0,1]. A rigorous mathematical treatment of the 
problem (3) for general classes of the weights fi was given in 115] . Some nonlinear equations 
with distributed order derivatives were studied in [I]. As above, this kind of relaxation models 
is connected with models of ultraslow diffusion [U [TUl [121 [HI [IZl [22]; these papers contain 
further references on related subjects. 

Let ux(t) be the solution of the problem (3) (the notation is changed slightly, compared to 
[H]). If /i(0) ^ 0, then 

MA(t)~C(logt)-\ t^oo. (5) 
If /i(a) ~ aa'^, a ^ (a > 0, > 0), then 

^ZA(^)~C(log^)-l-^ t^oo. (6) 

It was assumed everywhere in [2] that /i G C^[0, 1] and /i(l) 7^ 0; in fact, in the investigation 
of the asymptotic behavior of ux we use only that /i G Li(0, 1). Therefore the arguments in 
[13] cover the case where 

fi{a) ~ aa^'^, a ^ 0, 
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with a > 0, < u < 1, and yield the asymptotics 

^/A(t)~C^(log^)-^+^ t->oo. (7) 

In all the above cases, the exponential function e~^^, the Mittag-Leffler powerlike evolution 
Eai—Xt"), and all the logarithmic evolutions (5)-(7), the resulting functions are completely 
monotone, that is (— > 0, j = 0, 1,2,.. ., for all t. 

In this paper we look for other possible relaxation patterns corresponding to the weight 
functions tending to (at the origin) faster than the power function (Section 2) or to the 
definition of the distributed order derivative not by the formula (4) but by the expression 

1 

(D(,)ti)(t) = j (t)dp(«), (8) 



where p is a jump measure (Section 3). In these cases the solutions remain completely monotone 
while their asymptotic behavior can be quite diverse, from the iterated logarithmic decay to a 
function decaying faster than any power of logarithm but slower than any power function. 



2 THE MAIN CONSTRUCTIONS 

Suppose that p is a positive finite measure on [0, 1] not concentrated at 0, and consider the 
distributed order derivative (8). Substituting (2) into (8) we find that 

t 

d 



it) = - I k{t - T)u{r) dr - k{t)uiO) (9) 





where 

1 



Hs) = I ———dp{a), s > 0. (10) 







The right-hand side of (9) makes sense for a continuous function u, for which the derivative 

d * 

— j k{t — t)u{t) dr exists. 



dt 



It is clear from (10) that k G L'|°'^(0, oo), and k is decreasing. Therefore the function k 
possesses the Laplace transform 

oo 1 

]C{p)= jk{s)e-P'ds = jp'^^^dpia), Rep>0. 



The holomorphic function ]C{p) can be extended analytically onto the whole complex plane cut 
along the half-axis M_ = {Imp = 0,Rep < 0}. Obviously, IC{p) — > 0, as \p\ — > cxo (a precise 
asymptotics is found for some cases in pj^ and Section 3 below). 
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Note that if we consider the moments of the measure p, 

1 



and introduce their generating function 

1 

OO ^ n 

M{z)=Y,mJ-^ = / e-'dp{x), 

n=0 { 

then /C(p) = p~^M{logp). The moment generating functions were studied by a number of 
authors (for example, [5| 120]). 

Considering the relaxation equation 

(t) = -Xu{t), t > 0, (11) 

with A > 0, we apply formally the Laplace transform (which is justified post factum, using the 
smoothness properties and the asymptotic behavior of the solution). 

For the Laplace transform ux{p) of the solution u\{t) of the equation (11) satisfying the 
initial condition ux{0) = 1 we get the expression 

Since pfCijp) — > oo, as p — oo, we have u\{p) ~ p^^, p = a + ir, a,r eM, \r\ oo. Therefore 
[6], -ua is the Laplace transform of some function ux{t), and for almost all t, 

7+joo 

d 1 f e'P^ Up) , , , 

7— ioo 

Note that for p E C \ M_ we have 

1 

Imp/C(p) = / IpI" sm{a argp) dp{a), 







so that lmp)C{p) = only for argp = 0. This means that p}C{p) + X 7^ 0, and the representation 
(13) is valid for an arbitrary 7 > 0. 

As in [14j, in the present more general situation we deform the contour of integration, and 
then differentiate under the integral, so that 

1 f ^pt ^(P) 



S- 
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where the contour 5"^^^ consists of the arc 

T^,u = {peC: bl = 7, I argpl < un}, ^ < a; < 1, 

and two rays 

^t'^ = G C : I argpl = ±ujit, \p\ > 7}. 

Just as in the case of a measure p with a smooth density considered in [ll] , it is easy to show 
that the function Ux belongs to C°°(0, 00) and is continuous at the origin; from its construction 
and the formula (9), it follows that the initial condition ux{0) = 1 is indeed satisfied. 

Let us consider first the case of a measure p with a continuous density, dp{a) = da. 
This case was investigated in [M], and it was proved that ux is completely monotone (various 
additional assumptions made in [H] and needed for other problems studied in that paper, were 
not actually used here). 

In this paper we consider the case of a different behavior of the density p near the origin, 
implying a different asymptotics of ux{t), as t ^ 00. 

Theorem 1. If p e C[0,1], and 

p{a) ~ aa'^e"", as a ^ 0, 
where a > 0, 7 > —1, (3 > 0, then 

MA(t) ~C(log^)-^-te-2^(l°s*)^ t^oo. (15) 



Proof. Let us write IC{p) as 

00 

)C(p) = p^^ / e'"^^ pAa) da, 2; = log-, 

J P 


where pi is the extension of p by zero onto M+, and use an asymptotic result for Laplace 
integrals from [19] (Theorem 13.1, Case 9). We get 

-1 



K,{p) ~ 2ap-^ ^-j K^+i{2,/Pz), p ^ +0, 

where, as before, z = log ^ (— > 00), Km is the McDonald function. It is well known that 
Km{t) ~ (f)^^^^"^^^e-*, t 00, so that 

}Cip)r^Cp-'L{-), p^+0, (16) 
P 



where 

L{s) = (logs)~i~3e 



2_3 _2v/;3(logs)3 
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It follows from (16) (or directly from the definition of K,{p)) that pK,{p) 0, as p — ^ +0. 
Therefore by (12), 

u\{p) ~ \^^IC{p), as p — > +0. 

Since we already know that the function u\ is monotone, we may apply the Karamata-Feller 
Tauberian theorem (see Chapter XIII in [9]) which implies the desired asymptotics of ux{t), 
t — s> oo. ■ 

In the case under consideration, the function ux{t) decreases at infinity slower than any 
negative power of t, but faster than any negative power of logt. It is also seen from (15) that 
the decrease is accelerated if > becomes bigger, so that the less the weight function ^ is 
near 0, the faster is the relaxation for large times. 



3 THE STEP STIELTJES WEIGHT 

Let us consider the case where the integral in (8) is a Stieltjes integral corresponding to a right 
continuous non-decreasing step function p{a). In order to investigate a sufficiently general 
situation, we assume that the function p has two sequences of jump points, /?„ and z/„, n = 
0, 1,2,.. ., where /?„ — >• 0, — >• 1, Pq = uq E (0, 1). We may assume that the sequence {/?„} is 
strictly decreasing while {i^n} is strictly increasing. 

Denote Ap{t) = p{t) — p{t — 0), = ApiPn) {n > 0), r]n = Ap(z/„) {n > 1); we have 
Cn, Vn > for all n. It will be convenient to assume that /5o < and to write rjQ = 0. Since p 
is a finite measure, we have also 



J^^n<00, ^r]n <00. 



(17) 



n=0 



n=0 



By (10), 



n=0 



n=l 



S > 0, 



SO that 



n=0 n=l 

As before, we denote by u\{t) the solution of the relaxation equation (11) with the initial 
condition 'Ua(O) = 1. The symbol f g will, as usual, mean that / = 0{g) and g = 0{f). 



Theorem 2. (i) The function u\ is completely monotone. 



[II) Ux[X) 



n=0 



Vn 



[V{2 - I3n) r(2 - 

:iii) // Y: ?n (fog fog ^) <^ (h> O;, then 

1 



-X 



X ^ oo: 



(19) 



U\[X) 



O 



(log logx)^ 



X 



oo. 



(20) 
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(iv) // E ^nPn' <oo (b> 0), then 

n=0 



ux{x) = O 



X ^ oo. 



Proof. Using the representation (14) and following [13] we can write 



pfC{p) + A TT 



(21) 



— Ini 

TT J r (re*'^'^/C(re*'^'^) + A) 



def 



dr = Ji + J2 - Js- 



Let us substantiate passing to the limit as 7 ^ 0. 
As p — > 0, plC{p) — > 0, so that 



/C(p) 



p/C(p) + A 



n=0 



whence 



n=0 



as 7 — > 0. It was shown in [Til that 



J2 



TT 



— I s * sinfstancjyr) (is, 



as 7 — s> 0. 

The integral J3 is the sum of 



h 



Im 



Re 



1 



r / \re*^''/C(re*^'') + A 



and 



We have 



Im 



Re 



Im 



^-L^trcosun sin(tr siu wvt) 



dr. 



and this expression has a finite limit, as r — > 0. Since also pJC{p) — >■ 0, as p — > 0, we see that 
we may pass to the limit in Ji, as 7 ^ 0. 
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Let ^ 

$(r, a;) = Im ■ 



Substituting (18) and denoting 



,n=0 



f oo 

+ \ [^y^ sin(a;7r/3„) + r^^r^^" sin(a;7rj/„)] 



n=0 



we find that 



£ [Cn^^" sin(a;7r/3„) + T^^r''" sin(a;7ri/„)] 
$(r,a;) = -"^ 



G{r,u) 

The denominator tends to A^, as r — > 0. Noting that 



Rc =r"'e*""°^'"^cos(irsina;7r) 

and using (17) we find that the integrand in I2 belongs to i^i(0, 00). 
Passing to the hmit 7 — > we obtain the representation 

00 00 

u^{t) ^ -- j s~^e~*sin(stana;7r) ds j r~-^e*''''°^'^''sin(trsina;7r)*(r,a;) dr 



00 

j r"^e*''^°^^'^cos(trsinu;7r)$(r,a;)cZr 


where 

00 

E [inr'^'" cos(a;7r/?„) + Ty^r''" cos(u;7rt'„)] + A 

*fr a;) = — 

^'^ G{r,uj) 

Next, let us pass to the hmit, as a; — > 1. It follows from the Lebesgue theorem that the first 
two integrals tend to zero, and we get 



00 

^ 7 E [^nT^^ sin(7r/?„) + 77„r^" sin(7ri/„)] 



that is, up to a positive factor, u\ is the Laplace transform of a locally integrable non-negative 
function bounded at infinity. Therefore ux is completely monotone. 
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(ii) It will be convenient to turn to the Laplace-Stieltjes transform instead of the Laplace 
transform. Set 

X X 

= j k{s)ds, vx{x) = J ux{s)ds, 



so that the Laplace-Stieltjes transforms are as follows: 

oo 

= J e^^"" dx{x) = /C(p), vx{p) = u\{p). 



By (12), vx(p) = l(p)>i(p) where l(p) = is a slowly varying function near the origin. 

plC{p) + X 

We have 



Ef(^^- + Efp^^'-. -0. ,22) 



x[x) = 

n=0 

The function x is monotone increasing. Denote 

x*(C) = limsup^, C>1. 

x—too >C\X ) 

Since (3n,Vn G (0, 1), we find that 

x(Cx)<Cx(x), (23) 

so that x*(C) ^ C- Thus, x is an 0-regularly varying function (see [2|, especially Corollary 
2.0.6). 

On the other hand, it is seen from (22) that x*(C) > 1 for ^ > 1, so that x*(+l) = 1. Thus 
we are within the conditions of the "Ratio Tauberian Theorem" ([2j, Theorem 2.10.1), which 
yields the relation v\{x) ~ /(i)x(x), x — >• oo, so that 

v\{x) ~ Cx(x), X — oo. (24) 

In order to pass from (24) to (19), we have to check further Tauberian conditions. Since u\ 
is completely monotone, it is, in particular, non-increasing, thus belonging to the class BI (see 
Section 2.2 in [2] for the definitions of this class and the class PI used below). It follows from 
(23) that also x G BI. 

Next, 

oo J. 



n.=0 



Since > for C > 1, we have 



El rfl—B,^ _|_ Vri rfl — Vr, \ 

n=Q ^ ' 



-\-v\ ) 1 "■=2 



oo 

r(2-i/„)'^ 



OO 

I r(2-/3„)'^ ^ r(2-i/„)-^ 



n=0 
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where 



oo 



n=0 ^ ' 

as x — > oo. Thus, x G PI. 

Now we are within the conditions of the O- version of Monotone Density Theorem ([2], 
Proposition 2.10.3), which imphes the required asymptotic relation (19). 

(iii) — (iv) Let us prove (20). The proof of (21) is similar and simpler, and we leave it to a 
reader. It is obviously sufficient to deal with the first summand in each element of the series in 
(19). 

Let us consider the function 

Lp{x) = x~"'{\og\ogxY ^ X > e, (25) 

where a,b > 0. We have 

ip'{x) = x~"'{loglogxY~^-ip{x) 

where ip(x) = — aloglogx. It is easy to check that ip decreases on [e,oo), '0(e) = b, 
il){x) — ^ — cxD, as X — > oo. The maximal value of the function Lp is attained at a single point Xq 
where iP{xq) = 0, that is 

= a log log Xq. 

log Xq 

Denote yo = logxo- Then — = alogyo) so that yologyo = —. We will in fact need an 

Vo a 
asymptotic behavior of yo as a function of a, as a — ^> 0. It is known (see Section 1.5.2 in [8]) 

that 

ft 6 /loglogf 

yo = log - - log log - + O — ^ 

a « V 

Therefore ^ ^ 

Ci- flog -] <xo< C2- (log -] (26) 
a \ a J a \ a J 

where the constants do not depend on a, h. 
We have tp{x) < f{xo), so that, by (25), 



(log logx)'' 

We need an estimate of (^{xq) making explicit its dependence on a. By (26), 



a / \ a 



- ] = e-^i^g^ ^ 1 and (log = e"'°si°f5^ ^ 1, as a ^ 0. Next, Xq < — , so that 

a J \ aj ^ 

log log Xq < log(log C4 + log -) < log(C5 log -) < C& log log - 

a a a 
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whence f{xQ) < C7(loglog ^)^. 
As a result, 



(log log x)**' 



and we have proved (20). 
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